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Abstract. I argue that the rotation of white dwarfs is not 
a remnant of the angular momentum of their main sequence 
progenitors but a result of the mass loss process on the AGB. 
Weak magnetic fields, if present in stellar interiors, are likely 
to maintain approximately uniform rotation in stars, both on 
the main sequence and on the giant branches. The nearly uni- 
form rotation of the core of the Sun is evidence for the ex- 
istence of such fields. Exactly axisymmetric mass loss on the 
AGB from uniformly rotating stars would lead lead to white 
dwarfs with very long rotation periods (> 10 yr). Small ran- 
dom non-axisymmetries (~ 10 -3 ) in the mass loss process, on 
the other hand, add sufficient angular momentum to explain 
the observed rotation periods around one day. The process il- 
lustrated with a computation of the probability distribution of 
the rotation periods under the combined influence of random 
forcing by weak nonaxisymmetries and angular momentum loss 
in the AGB superwind. Such asymmetries can in principle be 
observed by proper motion studies of the clumps in interfero- 
metric images of SiO maser emission. 



Key words: stars: rotation - magnetic fields 
- stars: evolution - stars: mass loss 



white dwarfs 



1. Introduction 

White dwarfs are observed to rotate with typical periods of 
a day. The main sequence progenitors of these stars are also 
rotating, and it is generally assumed that the rotation of white 
dwarfs is a remnant of this main sequence rotation. Arguments 
involving conservation of angular momentum can be used to 
make this plausible (e.g. Perinotto 1990, Pijpers 1993). A prob- 
lem with this picture is, however, that progenitors of WD have 
gone through a giant stage in which at least the envelope ro- 
tated very slowly. Thus, it is necessary to assume that the cores 
of giants remain decoupled rotationally from their envelopes 
during the entire evolution from main sequence turnoff till the 
formation of the WD. Since little is known with certainty about 
the processes that might redistribute angular momentum in- 
side stars, this assumption can not be easily rejected. 

There are, however, observational and theoretical reasons 
to doubt this picture. A strong observational argument is the 
internal rotation of the Sun. The most recent helioseismologi- 
cal measurements (Elsworth et al. 1995, Kosovichev et al. 1997, 
Corbard et al. 1997) show that the rotation below the convec- 
tion zone is esentially uniform, with measured degrees of dif- 



ferential rotation well below the 30%. seen at the surface. The 
known hydrodynamic angular momentum transport processes, 
even with rather optimistic estimates of their efficiency, leave 
the Sun with a much too rapidly rotating interior (Spruit et 
al 1983). A new hydrodynamic mechanism recently studied in 
some detail is friction by internal gravity waves excited by the 
convection zone (Press 1981, Spruit 1987, Zahn 1990, Schatz- 
man 1993). Realistic calculations of this process appear to be 
difficult, but estimates indicate that it can be more effective 
than the other hydrodynamic processes (Zahn et al. 1997). 

Magnetic fields, on the other hand, have long been known 
to be very efficient at transporting angular momentum. The 
torques exerted by magnetic fields become significant already 
at very low field strengths. For the Sun, for example, a field 
of less than 10G can provide sufficient torque to maintain the 
observed uniformity of rotation. A number of mechanisms can 
provide such weak fields, for example a fossil field (remnant of 
the star formation process) or a dynamo- like process operating 
on (a small remnant of) the differential rotation of the core. 

In this paper, I develop the consequences of assuming that 
the cores of giants do, in fact, corotate approximately with 
their envelopes. After discussing the observational evidence on 
WD rotation rates I develop theoretical arguments for the ex- 
istence of effective coupling between the core and the envelope. 
This predicts very slowly rotating cores in the giant progenitor 
of a single WD. The rotation of single white dwarfs must then 
be explained by other processes. 

The same applies to neutron stars born in red giants. The 
observed pulsar rotation periods of the order of a second are 
much shorter (by a factor 10 3 or so) than expected if they 
formed in approximately uniformly rotating giants, and with 
our assumption of strong coupling of the core another mech- 
anism also has to be found to explain the rotation of pulsars. 
The processes differ somewhat for white dwarfs and neutron 
stars. The arguments for the neutron star case are developed 
in a companion paper (Spruit and Phinney, 1998). There, we 
show that the kicks with which neutron stars are born (as in- 
ferred from their transverse velocities) also impart angular mo- 
mentum at amounts sufficient to explain the rotation of most 
pulsars. 

To explain the typical rotation rates periods of single white 
dwarfs (of which only about 30 have measured rotation rates) , 
I show in Sect. [E| that small asymmetries in the mass loss pro- 
cesses during the last phases of evolution on the AGB are suffi- 
cient to explain the observed rotation rates. These asymmetries 
act as a random forcing through which angular momentum ac- 
cumulates in the envelope. A balance results between this ran- 



2 



dom forcing and the loss of angular momentum by the wind. 
The evolution of the angular momentum as the mass loss pro- 
ceeds turns out to be mathematically the same as that of the 
velocity of a particle experiencing Brownian motion in a gas, 
and can be described by a Fokker-Planck equation. Solutions 
of this equation (Sect. |5.l|) show that probability distribution 
of the angular momentum is close to a Maxwellian. The mean 
angular momentum decreases as the square root of the enve- 
lope mass remaining. Current observational evidence relating 
to the asymmetries needed in this picture is discussed in Sect. 

I 

1.1. Rotation speeds of AGB cores 

Starting with a rapidly rotating main sequence star, and as- 
suming uniform rotation during the expansion to the giant 
stage, we can estimate the rotation periods to be expected for 
white dwarfs evolving from single stars. An early type main se- 
quence (MS) star, rotating near its maximum speed (of the or- 
der 400 km/s), and expanding without angular momentum loss 
onto the asymptotic giant branch (AGB), has a rotation period 
P G = 2tv/{GM/R g ) 1 ^(R g /R ms ) 1 ^ ~ lOyr for R G ~ 1AU 
(except for a modest difference in gyration radius neglected 
here). Most early type MS stars rotate significantly slower, so 
that periods of the order 30-100 yr would be expected for the 
AGB descendants of early type stars. Some of the observed 
white dwarfs must have descended from solar type stars (F- 
G), which have periods of the order 30d at the end of their 
main sequence life. The AGB progenitors of these WD would 
rotate 100 times slower, with periods of the order of a thousand 
years. 

If the small amount of envelope mass is ignored which set- 
tles back onto the core during post- AGB evolution (more about 
this in Sect. [Oh , these rotation periods would also be inher- 
ited by the WDs formed. While there are a few magnetic white 
dwarfs with inferred rotation periods of at least a century, 
most WDs for which periods are known rotate much faster. 
We evidently need another mechanism to explain the rotation 
of typical single WDs. Before entering the discussion of possi- 
ble mechanisms, I briefly review the observational evidence on 
WD rotation. 

2. White dwarf rotation periods 

Three methods exist for measuring the rotation rates of white 
dwarfs. The largest number of determinations comes from spec- 
tral and/or polarization variations in magnetic WDs. Figure |l] 
shows these rotation periods, as compiled by Schmidt and Nor- 
sworthy (1991, see also Schmidt and Smith, 1995). Added to 
this compilation was G 158-45 (Putney, 1996) with a period of 
4.44 hr. A second group of determinations comes from ZZ Ceti 
and other oscillating stars. In a number of these, sufficiently 
detailed observations exist to identify the oscillation modes, al- 
lowing determination the period splittings due to rotation. The 
rotation periods for 7 oscillating white dwarfs collected from 
the literature (Table 1) are also shown in Fig ^. The spike in 
the figure at 3 10 4 d represents the (magnetic) stars whose pe- 
riods are inferred to be longer than a century, on the basis of 
the absence of variations in the polarization on time scales of 
decades. These stars were put at their approximate lower limits 
of 100 yrs. 

The widths of the narrow NLTE line cores have been used 
to set limits on rotation velocities of stars for which the mag- 
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Fig. 1. Rotation periods of isolated white dwarfs. Dark: oscillat- 
ing WD (asteroseismological periods), light: magnetic WD. From 
Schmidt and Norsworthy 1991, Putney 1996, and refs in Table 1. 
The peak at the right represents the lower limits for the 5 stars 
whose rotation period is larger than about a century. 



netic and seismological methods can not be used (Wesemael 
et al. 1980, Koester and Herrero 1988, Koester et al., 1998, in 
prep.). The detection limit, apparently around 20 km/s, is not 
sensitive enough to determine the rotation of stars in the ~ Id 
main peak in Fig. 1, but may be useful in setting limits on the 
number of rapidly rotating (P < lhr) stars. Reid (1996, his 
Sect. 3.1) and Heber et al. (1997) infer upper limits from 8 to 
40 km/s from Keck spectra of some 25 single white dwarfs. 

The distribution of periods in the main hump around Id 
looks the same for the magnetic and the oscillating stars, given 
the limited statistics. Very long periods are absent from the 
sample of seismologically determined periods, but this may be 
due to observational limitations. No stars have had their os- 
cillations followed long enough to detect period splittings of 
a decade. One concludes that with the (limited) data avail- 
able, there does not seem to be a significant difference in the 
distribution of rotation rates of magnetic (B > 10 5 G) and 
nonmagnetic (B <, 10 4 G) stars. There may, however, be other 
differences between the magnetic and nonmagnetic WD, apart 
from the field strength. Sion et al. (1988) and Liebert (1995) 
for example, argue that the magnetic stars are more massive 
than the nonmagnetic ones, and derive from more massive pro- 
genitors. 
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ref 


PG 1159-035 


1.38d 


Winget et al. 1991 


PG 2131+066 


5.07hr 


Kawaler et al. 1995 


PG 0122+200 


1.61d 


O'Brien et al. 1996 


GD 358 


1.38d 


Winget et al. 1994 


GD 165 


4. 2d 


Bergeron et al. 1993 


GD 154 


2.3d 


Pfeiffer et al. 1995 


G226-29 


8.9h 


Kepler et al. 1995 


Table 1. Asteroseismologically determined WD rotation peri 
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3. Angular momentum transport processes 

3.1. internal gravity waves 

The kinetic energy of convective motions in the envelope ap- 
pears as a source of pressure fluctuations at the boundary be- 
tween the envelope and the core. These fluctuations propagate 
as internal gravity waves through the core. The waves carry 
angular momentum and their dissipation therefore transmits 
torques. Assume that prograde and retrograde waves are ex- 
cited with the same amplitude. This is a good approximation 
if the rotation is slow, such that the effect of Coriolis forces on 
the wave generation process is small. As the waves propagate 
into the core, they conserve their angular momentum (or wave 
action) until dissipation becomes important. If the dissipation 
of prograde and retrograde waves is the same, no net angular 
momentum transport takes place. Prograde (retrograde) waves 
propagating into a medium of increasing (decreasing) rotation 
speed, however, meet critical layers (where the rotation rate 
equals the horizontal component of the phase speed), and are 
much more effectively absorbed there (see the discussions in 
Goldreich and Nicholson 1989, Zahn et al. 1997). Due to this 
asymmetry, there is a net angular momentum transport which 
tends to reduce the differential rotation. In effect, the internal 
gravity wave field is a source of friction. Zahn et al. find that 
for the Sun, the time scale for synchronization between core 
and envelope due to this friction is of the order 10 7 yr. Since 
this is of the order of evolution time scale to the giant branch, 
the process could be significant in maintaining synchronization 
during core contraction. Detailed evolution calculations includ- 
ing this process by Talon and Zahn (1998), however, still yield 
too large internal rotation rates for the present Sun. 

3.2. magnetic torques 

Magnetic torques are transmitted by the stress component 
B r Bj,/(4:n). Approximating this as constant over a spherical 
surface, the angular momentum balance is 



IQ = -BrB^r 6 



(1) 



The synchronization time scale between a core of radius 
r = 10 10 iiiocm, mass M — tuMq and moment of inertia / = 
k 2 MR 2 rotating at a rate Q = 10~ 6 fi_6S~ 1 is then 



3 k 2 MQ 
2 BrBfyR 



3 10 8 (B r B^ 



'• 10 1 n_ 6 



yr. 



(2) 



where k is the gyration radius. At the Sun's current rota- 
tion rate f2_e = 3, the synchronization time scale is less than 
the age of the Sun if (B r B^) x ^ 2 > 3G. If on its way to the gi- 
ant branch the core of the Sun were to contract to a radius of 
10 9 cm while conserving angular momentum, it would rotate at 
rate fi_6 ~ 300. To maintain corotation on the 10 7 yr contrac- 
tion time scale, a field strength (BrB^) 1 ^ 2 = 300G is sufficient, 
at this rotation rate. 

If the azimuthal and radial field components are of similar 
magnitude, (|l|) can be written in terms of the magnetic energy 
E B = B 2 /(8ir)(4ir/3)r 3 , and rotational energy En = ^IQ 2 as 



(3) 



If the B r and B$ are not comparable, (g) is only the mini- 
mum magnetic energy required. Nevertheless, since rfi is typ- 
ically such a large number, (|^) shows that corotation can be 
maintained by a magnetic energy which is small fraction of the 
rotational energy of the star, for spindown on a time scale long 
compared with the rotation period. 

3.3. Winding up of field lines 

Are such magnetic field strengths plausible? If the fields of the 
magnetic A stars are fossil (which unfortunately is still un- 
clear), sufficiently strong fields might also exist in the cores 
of solar type stars. Even if the initial fields (on the ZAMS) 
are lower than these values, however, differential rotation will 
increase the field strength quickly to values that have an ef- 
fect on rotation. Whether initially present in the star or de- 
veloping later by core contraction, differential rotation winds 
up the field lines, increasing the field strength. This problem 
has been studied in various forms since the '50s. Winding up 
of an axially symmetric poloidal field into a predominantly 
azimuthal field by differential rotation produces an opposing 
torque that is linear in the number of differential turns made, 
as in a harmonic oscillator. The result is an oscillation of alter- 
nate winding up and unwinding at a period given by the Alfven 
travel time through the star (Mestel, 1953), where the Alfven 
speed is that of the poloidal field (which is unaffected by the 
winding-up). Since Alfven waves travel decoupled from each 
other, each on its own magnetic surface, the oscillation period 
is different on each magnetic surface. The oscillations on these 
surfaces therefore get out of phase after a few oscillations, and 
the length scale across the surfaces decrease as t~ . In a finite 
time, dissipative processes across the surfaces become impor- 
tant, and the oscillation damps out by phase mixing (Spruit 
1987, Charbonneau and MacGregor 1993, Sakurai et al. 1995). 
The net effect of the process is that the component of differ- 
ential rotation along the field lines is damped out, on a finite 
time scale, and this can happen with an initial field that is 
much weaker than estimate 



3.4- Magnetic shear instability 

Another possibility is that a turbulent field is generated by the 
same magnetically mediated shear instability that has been 
shown to operate effectively in accretion disks (Hawley et al. 
1995, Matsumoto and Tajima 1995, Brandenburg et al. 1995). 
The conditions for magnetic shear instability to exist in a star 
have already been studied in detail by Acheson (1978, 1979) 
though the proper interpretation of this instability (Balbus and 
Hawley, 1992) was not clear at the time (see, however, Fricke, 
1969). In the context of stellar interiors, it has been studied 
again recently by Kato (1992), Balbus and Hawley (1994) and 
Urpin (1996). Wherever this instability exists it will lead to 
very rapid growth (on the differential rotation time scale) of 
a turbulent magnetic field, which then acts on the differential 
rotation like an effective viscosity. 

Acheson's (1978) analysis of the instability conditions in- 
cludes (unlike the more recent works) the effects of thermal 
and magnetic diffusion and of viscosity. The inclusion of ther- 
mal diffusion is especially important since it makes the insta- 
bility appear under much wider conditions. This is seen from 
Acheson's condition (7.27, a special case of his more general 
condition), which is equivalent to 
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where 



q = dlnfi/dlnr, F = dlnB /dlnr, 



(for instability), (4) 



(5) 



N is the buoyance frequency, 7/ and k the magnetic and 
thermal diffusivities, and H the pressure scale height. This 
condition holds for low viscosity {u/rj <JC 1), for an azimuthal 
field B$ at the equator of the star. The first term on the left 
hand side represents the magnetic shear instability, the second 
term Parker instability (magnetic buoyancy instability). For 
weak fields, this second term is negligible. The right hand side 
shows the stabilizing effect of the stratification, which, how- 
ever, is partially undone by thermal diffusion (for adiabatic 
perturbations, the factor tj/k would be replaced by unity). 
Since photons diffuse so much more effective than the magnetic 
field, the instability is present much more widely than in an 
adiabatic treatment. The instability, however, is able to grow 
only on length scales sufficiently small that thermal diffusion is 
important. This somewhat limits its effectiveness, and it may 
be that the effective viscosity it produces is not much larger 
than the viscosity produced by hydrodynamic shear instabili- 
ties (Zahn, 1974) under the same conditions. These questions 
could, in principle, be readily addressed by an appropriate nu- 
merical simulation. 

Because magnetic fields are so effective at transmitting 
torques, already at low field strengths, differential rotation can 
survive over a large number of rotations only in regions where 
the radial field component is very small. In order to allow the 
core in a giant to rotate substantially faster than its envelope, 
one must find a reason why it could have been so accurately 
'shielded' magnetically, over the entire life of the star on the 
giant branch. 

While the arguments given here do not constitute a proof, 
I feel they are sufficiently strong that approximately uniform 
rotation (modulo a factor of a few) is a reasonable hypothesis, 
compared with the traditional assumption in which the core of 
a giant rotates 10 4 -10 5 times faster than its envelope. 

4. Mass and angular momentum loss on the AGB 

A large fraction of the star's mass is lost in the last phases of 
evolution on the AGB (e.g. Habing, 1990). Most is ejected in 
the form of a superwind (~ lO _4 M0yr _1 ) lasting on the order 
of 10 4 yr (e.g. Vassiliadis and Wood, 1993). It is believed to 
be driven by pulsational instability and radiation pressure on 
dust (Fleischer at al. 1992, Sedlmayr and Carsten 1995, Hofner 
and Dorfi, 1997), or possibly by sound waves (e.g. Pijpers and 
Hearn, 1989). The mass loss is probably not steady because the 
stellar pulsation is an important part of the driving. Also, dust 
formation in the expanding flow is an unstable process (Hofner 
and Dorfi, 1997). Thus the envelope is probably ejected in the 
form of a (large) number of light shells. The mass loss is also 
believed to be modulated on longer time scales by the thermal 
pulses of the AGB star. 

A small fraction of the star's envelope (on the order 
1CP 4 M0) settles back onto the core after the superwind ceases. 
Most of the angular momentum is lost together with the mass 
of the envelope, but because of the large size of the envelope, 
even the small amount of mass remaining might conceivably 



contain enough angular momentum to form a significantly ro- 
tating white dwarf. Thus we need to look in some detail at the 
angular momentum balance of the mass losing AGB envelope. 
First, I show that if the superwind is axially symmetric and 
has the specific angular momentum of the stellar photosphere 
from which it is ejected, the angular momentum remaining af- 
ter envelope ejection is far too small to produce a significantly 
rotating white dwarf. 

4-1. Axially symmetric mass loss 

If the mass is ejected from the stellar photosphere in axisym- 
metric fashion, taking with it the angular momentum it had in 
the photosphere, the net angular momentum loss by the wind 



(6) 



where M is the mass loss rate, R,, fl, the photospheric 
radius and rotation rate of the envelope. The factor 2/3 is due 
to the variation of specific angular momentum over the surface. 
Since the envelope is convective, it is a good approximation to 
assume that it rotates uniformly. Because of the very large 
radius of the envelope, the core contributes very little to the 
star's moment of inertia, even if the envelope mass is quite 
small. By angular momentum conservation the star's angular 
momentum J» varies as 

j, = lM,R 2 A, (7) 

where M, is the star's mass. With uniform rotation, J* = 
k 2 M*QR 2 , where k is the radius of gyration, hence 



2M„ 
3fc 2 M, ' 



(8) 



In stars, k < 0.4, so that the angular momentum of the 
star decreases more rapidly than its mass. This is because the 
specific angular momentum of the mass leaving the star is 
higher than the average specific angular momentum of the star 
(by a factor 2/3fe 2 ). Since the envelope mass varies strongly, 
the gyration radius can not be taken as constant. The total 
moment of inertia of the star can be written as the sum of core 
and envelope contributions: 



J* = k 2 M, 



*R 2 = j P vj 2 A i r+ j 

J core J c] 



Ic+h 



i envelope 

where w is the distance to to rotation axis, and 



/e = k 2 McRc, h = k 2 MeRl 



(9) 



(10) 



If the envelope contains most of the stellar mass, k e is ap- 
proximately that of a polytrope of index 1.5, k 2 « 0.2. For the 
estimates below I assume this value. For a degenerate core of 
mass ~ 0.6, k 2 is of the order 0.19. The gyration radius of the 
star as a whole is then 



'< = I/(M,Ri) = k\ 



Me^£ ,2 Me 

M,R 2+ ° M, 



(11) 



For R* ~ 10 13 , R c ~ 10 9 , the first term is negligible for 
envelope masses larger than ~ 10~ 8 M Q , so that 
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With ft) this yields 



ill M * 
J./J.= m 



(Me J> 10~ 8 ), 



where 
2 



3.3. 
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(14) 



Since the core mass is essentially constant during the mass 
loss, we have M w M c . Eq. @ can be integrated to yield 



Jo VM e0 J ' 



where Jo and M e o are the initial angular momentum and 
envelope mass. The steep dependence on M e implies that only 
a small fraction of Jo is retained. An upper limit on the final 
rotation rate is obtained by assuming the AGB star to rotate 
critically, Q* = {GM/Rl) 1/2 . The rotation rate of the post- 
AGB core then becomes 



n k 2 M* Rl /GM,\^ 2 /M ef 



(12) 5. Slightly nonaxisymmetric mass loss 

The angular momentum evolution of the star is altered dra- 
matically if even a small amount of non-axisymmetry is al- 
lowed in the mass ejection process. If a shell is ejected as- 
pherically, it generally carries a net momentum, and the direc- 
tion of this momentum vector in general need not pass exactly 
through the center of mass. It is conceivable, for example, that 
the dust-formation instability found in spherically symmetric 
numerical simulations of the ejection process actually is non- 
axisymmetric, so that the forces exerted are not evenly dis- 
tributed over the surface. In this way, the ejection process adds 
a small amount of angular momentum ('kick') to the star. 

Suppose now that a large number of shells are ejected, 
adding small amounts of angular momentum in random direc- 
tions. Since the simulations indicate that the shells are ejected 
with periods St of the order of the oscillation period of the star 
(on the order of a year) while the duration of the superwind 
phase is of the order of 10 4 yr, there are on the order of 10 4 
kicks, each associated with an ejected mass SM on the order of 
1O -4 M0. The maximum amount of angular momentum such 
a kick can impart is 5MR 2 v e , where v e is the ejection velocity, 
observed to be in the range 5-50 km/s. This maximum applies 
when the mass is ejected tangentially to the surface of the star. 
This is of course quite unrealistic, and one expects the angular 
momentum imparted to be only a small fraction of this: 



(15) 



(16) 



where indices f and * denote the post-AGB core and the 
AGB star, respectively. If at the end of the superwind phase 
an envelope mass of not more than 10 -3 M o is left, we get a 
final rotation period of at least a year. 

The effect depends rather critically on the index m in dli|). 
If the wind corotates with the star out to some radius R > 
R t , for example because of an atmospheric magnetic field, the 
specific angular momentum in the wind is increased by the 
factor / = (R/R„) 2 , and the index m would become 



2/ 
3h c 



(17) 



Magnetic fields are known to exist in Mira envelopes from 
the circular polarization of the SiO masers (Barvainis et al. 
1987, Kemball and Diamond 1997). The values of the field 
derived are uncertain since they depend on the degree of satu- 
ration of the masers (Nedoluha and Watson 1994). A strength 
of a few tenths of a Gauss, however, would already cause sig- 
nificant additional angular momentum loss by the wind. 

The conclusion is that even a maximally rotating AGB star, 
with its huge amount of angular momentum, will produce only 
a nearly non-rotating white dwarf if mass loss is axisymmet- 
ric. This results from the fact that almost all the envelope is 
lost, combined with the higher than average specific angular 
momentum taken away by the mass lost. Physically, as mass 
is lost from the photosphere, the envelope expands, causing 
spindown by angular momentum conservation. 

Let me summarize the assumptions made in arriving at 
this, perhaps surprising, conclusion. The first is that core of 
the AGB star corotates approximately with the envelope when 
the phase of rapid mass loss sets in. The others are the rather 
minimal assumptions that the (convective) envelope rotates 
approximately uniformly, and that the mass lost in the wind 
carries at least the specific angular momentum of the photo- 
sphere of the star. 



8 J = eSMR 2 v c , 



(18) 



where e is a small number. In the following, I estimate how 
large this number must be to explain the observed rotation 
periods. 

The evolution of the star's angular momentum vector is 
obtained by adding the forcing by kicks to (|l3|): 



Me 



■A(t), 



(19) 



where A is a random fluctuating vector with time step St 
and amplitude SJ. This equation (Langevin's equation) is the 
same as that governing the Brownian motion of particles in 
a gas. Following the standard treatment in statistical physics 
(e.g. Becker, 1978) we can take the continuum limit, in which 
the time step is infinitesimal, and derive a Fokker-Planck equa- 
tion for the probability distribution /(J, t) of obtaining an an- 
gular momentum J after time t. Leaving out this derivation, 
the result is 



df/m = i3Vj-(Jf + DVjf), 

where 
p — m- 



Me ' 



(20) 



(21) 



is the 'braking rate', Vj is the gradient in J-space, and D 
the diffusion coefficient in J-space 



D= USjf/St. 



(22) 



The main difference with respect to standard Brownian 
motion is that the coefficients j3, D in the present case are 
functions of time. 

If the kicks are random in direction, and the star initially 
non-rotating, the probability distribution / is isotropic in J- 
space, / = /(J). Writing 
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J 2 f, 



Eq. (|2C|) can then be written as 
d t F = dj\fiJF + D(djF - 2F/J)]. 



(23) 



(24) 



If (3 and D are constant, as they are in the case of Brownian 
motion, the asymptotic solution F^o for large t is that for which 
the bracket on the RHS vanishes. This yields 



J 2 exp( 



2D J ' 



(25) 



i.e. a Maxwellian distribution peaking at J = (2D/ /3) 1 / 2 . In 
our case, j3 varies significantly with time, because the envelope 
mass varies strongly. We should therefore do not expect the 
distribution function to be a Maxwellian. 

Before entering into more detailed calculations an estimate 
of the orders of magnitude to be expected can be made by mak- 
ing a quasi-stationary approximation to (^) . The distribution 
F is then given approximately by jj^). The typical angular 
momentum to be expected at the end of the mass loss, when 
the envelope mass left is M a <C M c q is then, with (pl|), (|2^): 

J » SJ(nk 2 M c /M c0 ) 1/2 , 



(26) 



where 



n = (Meo - M E )/{M c St) (27) 
is the number of kicks experienced. With (|l^) : 



J « eM, 



,oR*v e ( 



V rnnMco 



(28) 



where m = 2/(3fc 2 ) for angular momentum loss at the pho- 
tospheric value (/ = 1 in ^) . The expected rotation rate of the 
white dwarf, with gyration radius fc w , mass M w and radius _R W 
is then Q w = J / {k% M w i? 2 ,) . Assuming a final envelope mass 
of 1O _4 M0, initial envelope mass of 2Mq, and M w = 0.6Mq, 
this yields 



e5n 



-1/2 
7 c 



(29) 



If shells are ejected every two years or so, we have n w 10 4 . 
A rotation period of Id is then obtained for e « 10 -3 . 

As long as it is not known how the relevant details of the 
ejection process take place, it is hard to argue whether an asym- 
metry of the order 10 -3 is realistic or not, but a number as 
small as this would not seem too demanding. The reason why 
such small asymmetries are sufficient, even when their effect is 
further reduced by random superposition (the factor n" 1 / 2 in 
|28| ), is the very large lever arm on which the kicks act. A star 
on the AGB is so large compared with the final white dwarf 
that a very precisely axisymmetric mass loss would be needed 
to avoid introducing the small amount of angular momentum 
that is sufficient to produce white dwarfs with periods of a day. 

5.1. Distribution of rotation rates resulting from random kicks 

The coefficient /3 varies by a factor 10 4 as the envelope mass 
is reduced from its initial value to a representative post-AGB 
value of the order 1O _4 M0. To take this into account, I solve 
Eq. (^) numerically. I use a second order, implicit time step 
and centered differences in the J-coordinate (Crank-Nicholson 
scheme) . 

As angular momentum coordinate I use the dimensionless 
variable j, defined by 



-0.5 




-3 -2 

log(M E /M E0 ) 



Fig. 2. The mean angular momentum j m as a function of the re- 
maining envelope mass M e , for angular momentum induced by ran- 
dom non-axisymmetries in the superwind mass loss. Solid: solution 
of Eq. (33). Dotted: stationary approximation, Eq. (34). Evolution 
is from right to left 



J = jeMeoR*v e (mn) 1 ' 2 . 
Let 

g(3)=fm=dN/dlnj 



(30) 



(31) 



be the probability distribution per unit of In j, and for time 
coordinate use 



r = -ln(M e /M e0 ). 

Then Eq. (|5j) can be written as 
d T g = jdj [mg + e~ T j 2 dj(g/j 3 )]. 



(32) 



(33) 



The integration is from r = to Tf = — ln(M e //M e o), 
where M e / is the envelope mass at which the mass loss ends. 
Apart from M / the only parameter in the problem is the an- 
gular momentum loss index m = 2//(3fc 2 ) (cf. Eq. ^). The 
evolution for m = 3.3 is given in Fig. H, which shows the mean 
jm of the probability distribution g(j, M e ) as a function of the 
remaining envelope mass. 

If the evolution is sufficiently 'slow', one expects the solu- 
tion to be close to the Maxwellian stationary solution, obtained 
by setting the square bracket in (^) equal to zero. This station- 
ary distribution has mean dimensionless angular momentum 



js = 2 V / 2MM C /M c0 ) 1/2 , 



(34) 



and is shown for comparison in Fig. ^. The stationary ap- 
proximation actually turns out to be quite good, except in the 
initial phase of the evolution. The white dwarf rotation rate 
corresponding to Q3-4J) is given by: 



- M c0 R, v c Mg 
fc|M w i? w Rw mnM e0 



\l/2 



(35) 



Thus, the predicted white dwarf rotation rate decreases as 
the square root of the mass remaining in the envelope at the 
time when mass loss ceases. 
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Fig. 3. Predicted distribution of rotation periods (solid) for asym- 
metry parameter e = 10 — 3 , and a log- normal spread in final envelope 
mass from 210~ 5 to 51O- 4 M . Histo gram: observed distribution 
from Fig. 1 



5.2. Comparison with observed distribution 



By adjusting either the asymmetry parameter e or the final 
envelope mass M e /, the maximum of the distribution can be 
made to agree with the observations. This distribution is close 
to a Maxwellian, and its width is too narrow compared with the 
observations, which spread by a factor 20 or so. The factors in- 
fluencing the mean rotation rate ^Ej) most are the asymmetry 
parameter e and the remaining envelope mass M c q. Both might 
depend on systematic factors like the initial stellar mass. A ran- 
dom variable could be the phase in the thermal pulse cycle at 
which the superwind takes place, which is known to have an ef- 
fect on the post-AGB evolution (Schonberner 1990, Vassiliadis 
and Wood 1993). Lacking a sufficiently detailed theory for the 
superwind, it is hard to guess how the asymmetry parameter 
might depend on such variables. Values of the remaining enve- 
lope mass, on the other hand, have been inferred for oscillating 
WD and post-AGB stars by comparisons with theoretical mod- 
els. Clemens (1994) finds a hydrogen envelope masses of about 
10 _4 Mq. In the helium (DB) white dwarfs and their possible 
progenitors the PG 1159 stars, only a helium envelope (with 
inferred masses of the order 1O _3 M , cf. Dehner and Kawaler 
1995) is left. Blocker and Schonberner determine a hydrogen 
envelope mass of 3 1O" 4 M for FG Sge. It seems reasonable 
to assume that a certain spread in M e / is present. This could 
be due, for example, to random variations in moment at which 
pulsation ceases. To fit the observed distribution with such a 
spread, I assume a log-normal distribution of the parameter 
M e /, with peak at M cf = 1O~ 4 M and (1/e-) width from 
2 10 -5 to 5 10 -4 . The asymmetry parameter is assumed to be 
e = 10 -3 . The resulting period distribution is compared with 
the observations in Fig. H. The agreement with the observa- 
tions is not a test of the theory developed here (since both the 
width and the mean have been fitted), but comparison shows 
that a spread in envelope mass of a factor 5 on both sides of 
the mean is sufficient to explain the observed width. 



5.3. Binarity and common envelope 

Alternatives to the picture sketched may be envisaged, in which 
the rotating white dwarfs are in one way or another the result 
of binarity. The following is a brief discussion of such scenarios. 

The observed periods, O.ld < P < lOd are suggestive of 
the orbital periods of close binaries, as has been noted by 
Schmidt et al. (1986), who suggest the possibility that AM 
Her stars might be the progenitors of the rotating white dwarfs. 
Those must then have somehow lost their companions, perhaps 
through the mass transfer. Current understanding of the evo- 
lution of CVs does not favor this possibility, since it predicts 
that the secondaries can not transfer all of their mass within a 
Hubble time. This is because the angular momentum loss slows 
down dramatically once the secondary has been reduced to a 
small degenerate dwarf (e.g. Verbunt 1996, Kolb 1993). 

A second possibility that suggests itself is that of a binary 
companion absorbed in a common envelope (CE) process. Two 
different outcomes of such a CE are possible. One is that the 
envelope is ejected, by the orbital energy released, before the 
secondary has spiraled in completely. In the other, the sec- 
ondary spirals in completely and merges with the primary. 
The first case leaves a detached system (such as V471 Tau) 
which then evolves into a CV by magnetic braking. Theory 
and numerical simulations (for reviews see Taam 1995, Livio 
1996) predict that this case happens if the secondary is mas- 
sive enough and the density gradient in the inner parts of the 
giant are not too steep. If these conditions are not met, the sec- 
ondary is predicted to dissolve completely, transferring all its 
mass to the giant envelope. A significant fraction of common 
envelope systems may actually experience this fate. 

The high incidence of elongated or bipolar structures in 
planetary nebulae (PN) and objects believed to be in transit 
from an AGB star to a PN suggests that a large fraction of 
PN involve some form of common envelope evolution (e.g. Han 
et al. 1995). Detailed hydrodynamical simulations have been 
made to reproduce the morphology of these nebulae (Icke et 
al. 1992, Frank and Mellema 1994). The results show that in 
the initial phases of the radiation driven nebular expansion 
there must have been a thick disk-like structure inhibiting fast 
outflow in the plane of the disk, leaving a structure of two 
rapidly expanding lobes and a more slowly expanding ring. In 
the CE interpretation, the disk contains the mass ejected in 
the spiral-in process. 

If the secondary is small, the energy released as it spirals 
in is insufficient to eject the entire envelope of the primary. 
The net effect in this case is that both the mass and the angu- 
lar momentum of the secondary are added to the envelope of 
the primary. The envelope remaining on the primary after the 
CE would then contain a large amount of angular momentum, 
even if the companion absorbed is small. Would this suffice to 
produce a rotating white dwarf? If our basic assumption of ap- 
proximately uniform rotation is valid, the answer is negative. 
This follows from the analysis of Sect. Ej] , where I have shown 
that even a maximally rotating AGB star leaves a core rotating 
with a period of at least 10 years. 

This answer applies as long as there is still a significant 
amount of mass left in the envelope after the common envelope 
process (O.1M0, say), and mass loss then continues like in 
normal AGB stars. If any significant amount of mass is left in 
the fo rm of a convective envelope after the CE, the results from 
Sect. 4.1 predict that the result will be a very slowly rotating 



white dwarf. 
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The consequence of the above is that a rapidly rotating 
white dwarf by CE evolution is to be expected only if the final 
dissolution of the companion coincides rather precisely with the 
ejection of the last bits of envelope. Barring possible surprises 
concerning late phases of CE evolution, the details of which 
are not well known, this situation would appear to be a rare 
coincidence. 



6. Discussion 

In the standard view, the rotating single WD derive from the 
rotating cores of giants, which somehow avoided spinning down 
in the slowly rotating convective envelope. I argue, instead, 
that rotating cores in giants are an unattractive idea, espe- 
cially if these cores are magnetic. Unless the magnetic WD 
acquired their fields after emerging from the envelope, the ob- 
served dipole moments are so large that a strong interaction 
with the slowly rotating convective envelope would be very 
hard to avoid. 

I recall the classical demonstration (e.g. in Mestel 1953, 
1961) that rather weak magnetic fields (magnetic energy a 
small fraction of the rotational energy) can already transmit 
enough torques to maintain corotation between core and enve- 
lope. Such a weak field could be inherited from the star forma- 
tion process. In order to prevent these torques from acting, any 
magnetic field in the core would have to be very weak or very 
accurately shielded from the convective envelope. In addition, 
a differentially rotating, initially nonmagnetic core is unstable 
to the growth of a small scale dynamo magnetic field, initiated 
by a magnetic shear instability (Balbus and Hawley 1992). The 
conditions for existence of this instability in stars were studied 
in detail already by Acheson (1978) who showed, in particular, 
that thermal diffusion allows it to operate under a much wider 
range of conditions than in the adiabatic case. 

The very weak differential rotation in the core of the Sun 
(e.g. Kosovichev et al.1997), for which no good explanation 
has been put forward except magnetic torques, is strong ev- 
idence for the operation of magnetic effects. While the argu- 
ments given here do not constitute a proof, I feel they are suf- 
ficiently compelling that approximately uniform rotation is a 
reasonable hypothesis, and is at least as plausible as the tradi- 
tional assumption, which implies a core rotating 10 4 -10 5 times 
faster than the envelope for the entire duration of the RGB 
and AGB. 

I have explored the consequences of the assumption of ap- 
proximately uniform rotation for AGBs stars in the process of 
shedding their envelopes. If this mass loss is strictly axisym- 
metric, the remaining core rotates very slowly (period more 
than 10 years). This is just the consequence of angular mo- 
mentum conservation: the wind takes away almost the entire 
envelope, but the specific angular momentum it carries away is 
that of the stellar photosphere, which is larger than the average 
specific angular momentum of the envelope. 

On the other hand, only small nonaxisymmetries in the 
mass loss process suffice to give the star enough 'kick' to ex- 
plain the angular momentum of single white dwarfs. Such kicks 
could be associated with mass loss events at the pulsation pe- 
riod of the star or dust-formation episodes in the atmosphere. 
I have illustrated this with a calculation of the evolution of the 
probability distribution of the star's angular momentum un- 
der the combined action of many small nonaxisymmetric kicks 



and the angular momentum loss in the wind. The degree of 
asymmetry required is found to be of the order 1CP 3 . 

Present theories for AGB mass loss are not detailed enough 
to calculate such asymmetries, but observational indications 
for asymmetries exist. Interferometric images of red super- 
giants (a Sco, q Ori and a Her: Tuthill, Haniff and Baldwin 
1997), speckle reconstructions (a Ori: Kluckers et al. 1997) 
and HST imaging (a Ori: Gilliland and Dupree, 1996) show 
pronounced 'hot spots' on their surfaces. Assuming that such 
nonuniform photospheric conditions persist during the super- 
wind phase, one would expect them to also affect the dust 
formation that is thought to be essential for the driving of the 
wind. The required asymmetry is obtained if a few (5 say) such 
spots are present, and the wind locally generated above these 
spots is slightly non-radial by a few tenths of a degree. That the 
mass flow is indeed asymmetric already close to the stellar pho- 
tosphere is shown by speckle imaging (IRC 10216: Osterbart 
et al. 1997), and especially by mm- wave interferometric images 
of the SiO maser emission. These show a highly clumpy and 
time dependent structure (Diamond et al. 1994, Humphreys 
et al. 1996, Pijpers et al. 1994). This maser emission occurs 
at a few stellar radii, which is also the region where the back- 
reaction of the wind on the star ('kick') takes place. Though 
the SiO maser emission is very sensitive to small changes in 
the local physical conditions, models of the emission (Lockett 
and Elitzur 1992, Bujarrabal 1994) should give estimates of 
the degree of inhomogeneity in the physical conditions in the 
wind. 

Measurement of deviations from radial flow in proper mo- 
tion studies of the masing clumps in the wind should enable 
direct determination of the asymmetries relevant for the kick 
process described in this paper. 

An issue mentioned here only briefly is the origin of the 5 or 
so very slowly rotating (P > lOOyr) magnetic white dwarfs. A 
possible explanation is angular momentum loss in a radiation 
driven, but magnetized, wind during post- AGB evolution. This 
possibility will be further explored elsewhere. 

The coupling between core and envelope proposed here 
would also imply that the cores of pre-supernovae on the giant 
branch are so slowly rotating that very slowly (P ~ lhr) rotat- 
ing neutron stars would result even if angular momentum were 
conserved during core collapse. While these would not show up 
as pulsars, one would have to argue that none of the observed 
pulsars were formed in red giants, which feels like an unattrac- 
tive idea. It turns out, however, that the kicks neutron stars 
receive at birth and which give them their high observed space 
motion, are strong enough to impart a significant rotation as 
well. This idea is developed further in a separate paper (Spruit 
and Phinney, 1998). 
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